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Abstract We study a specific correction to the Bethe logarithm in- 
duced by potentials which are proportional to a Dirac-<5 function in 
coordinate space ("local potentials"). Corrections of this type occur 
naturally in the calculation of various self-energy corrections to the 
energy of bound states. Examples include logarithmic higher-order 
binding corrections to the two-loop self-energy, vacuum-polarization 
induced corrections to the self-energy and radiative corrections in- 
duced by the finite size of the nucleus. We obtain results for excited 
S and P states and find that the dependence of the corrections on 
the principal quantum number is remarkable. For the ground state, 
we find a small modification as compared to previously reported re- 
sults. Our results are based on mathematical techniques for the treat- 
ment of quantum electrodynamic bound states discussed previously in 
[J. Phys. A 35, 1927 (2002), hep-ph/0111084]. 
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1 Introduction 



Ever since the quantization of the electromagnetic field was introduced by Dirac in [1], one 
of the main directions of research in the area of quantum electrodynamics has been the study 
of bound states and, notably, the corrections to the energy of these states as induced by the 
virtual quanta. A severe limitation to the accuracy of current theoretical predictions is given 
by bound-state two-loop self-energy effects whose evaluation has been historically problematic 
(see [2, 3] and references therein). In this Letter, we present complete results for the two-loop 
(2L) logarithmic self-energy correction of excited S states of order 

A< L> («S) = (i!)' (^HtZaJ-lBn^), (!) 

where Bqi(uS) is an re-dependent, dimensionless coefficient. Here, re is the principal quantum 
number, Z is the nuclear charge number, m is the electron mass, and a is the fine-structure 
constant (we work in natural units: fi = c = eq = 1). 

One of the quantities which enter quite universally into higher-order corrections to the one- 
and two-loop self-energy are those induced by effective potentials which are proportional to 
a Dirac delta-function in coordinate space. In order to calculate Bq\ for excited S states, an 
investigation of such corrections is necessary [3]. Our investigations are based on mathematical 
techniques which facilitate the treatment of bound states which rely on a separation of the 
virtual photons into hard (high-energy) and soft (low-energy) quanta; these have recently been 
generalized to two- loop effects [2, 4]. 



2 General Formulation 

Let us consider the Schrodinger Hamiltonian 

*-£ + V, (2) 

where V is the binding Coulomb potential (the energy eigenvalues are E = —(Za) 2 m/(2n 2 )). 
We assume a small perturbation of H proportional to a 5-function, 

For S states, this perturbation leads to an energy shift 

8E = (<j>\5V\<l>) , (4) 

where \<j>) is the electron wave function for which the nonrelativistic approximation may be used 
in the context of the evaluation of the radiative corrections as discussed in this Letter. The 
correction (4) is nonvanishing only for S states, where the correction to the raS-state energy 
amounts to (Za) 4 /n 3 . The perturbation (3) also induces a modification of the wave function 

where the prime denotes the reduced Green function (in the spectral decomposition, the reference 
state \(f>) is excluded). 
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We calculate the delta-like correction to the one-loop Bethe logarithm, 



2a 
3^ 



P 



P 



m E — (H + to) m 



-{Zaf ^F 6 (nl,e) 



(6) 



where nl denotes the bound-state quantum numbers. The notation (6) follows the usual spectro- 
scopic nomenclature, n is the principal quantum number, and / is the orbital angular momentum. 
The correction (6) is independent of the electron spin, and e denotes the quantity 

(7) 



(Za) 



(the notation has been introduced in Ref. [3]). The upper index (L) in Eq. (6) is assigned 
because AE^ represents the low-energy part of the correction (due to soft virtual photons), in 
the language of the formalism introduced in [2, 5, 6]. In Eq. (6), the symbol Sgy, inspired by the 
notation of [3], denotes the first-order perturbation received by the quantity in curly brackets 
through the replacements [see Eqs. (3), (4), (5)] 



H 



H + 5V 



(8) 



\<f>) -> \cb) + \5<P), (9) 

E -> E + 5E. (10) 

The quantity FQ(nl,e) in Eq. (6) is a dimensionless function which parametrizes an effect of the 
order of (Za) e . Here, we present results for nS states in the range n = 1, . . . , 8 and nP states 
(n = 2,...,8). 



N(nS) 




1 2 3 4 5 6 7 8 



Figure 1: The dependence of the correction iV(nS) on the principal quan- 
tum number is remarkable. For n = 8, the magnitude of the result is 
less than half iV(lS). The smooth curve is a fit based on a model based 
on a three-parameter fit of the form as + bs/n + cs/n 2 [see Eq. (14)]. 
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3 Correction to S states 



For S states, Fq as implicitly defined in Eq. (6) has the form [3] 

2{1 - In 2} 



F 6 (nS,e) = -| ln 2 (e)+ln(e) 



+ l (z + 4^4- ln(n)+ * (n)+c 



+ N(nS). (11) 



Here, the notation e is explained in Eq. (7), and iV(nS) is a nonlogarithmic term which has been 
known only for n = 1 (see Ref. [3]). The dependence on e cancels when the high-energy part is 
added to the above result [5, 6, 7]. The constant term iV(nS) has been shown to contribute to 
the two-loop self-energy coefficient Bqi(uS) [see [3, Eq. (50)]]. The complexity of the calculation 
increases with increasing principal quantum number, because of the more complex structure of 
the bound-state wave function and the necessity to subtract poles of the integrand corresponding 
to the decay into lower-energy states. For states with n = 8, we obtain intermediate results with 
198,000 terms, and use is made of computer algebra systems [8]. 

We obtain the following results for N(nS) is the range n = 1, . . . , 8: 



N(1S) 


= 17.855 672(1), 


N(2S) 


= 12.032 209(1) 


N{3S) 


= 10.449 810(1), 


N{AS) 


= 9.722413(1) , 


N{5S) 


= 9.304114(1), 


N(6S) = 


- 9.031832(1) , 


N{7S) 


= 8.840123(1), 


N(8S) = 


8.697639(1) . 



(12) 

A least-squares fit with a functional form 

N(nS) « a s + b s /n + c s /n 2 (13) 

yields the fit-parameter values 

a s = 7.78, 6 S = 3.13, c s = 6.93. (14) 

The data in Eq. (12) are well represented by this fit, as is evident from Fig. 1. The excellent 
agreement between the fit and the numerical data in Eq. (12) could suggest a finite limit 

lim N(nS) w 7.78. (15) 

n— too 

The functional form of the fit is inspired, in particular, by the term proportional to ln(e) in (11). 
This term can be expanded in a series involving inverse powers of n. Note that the logarithmic 
term ln(n), for large n, cancels against a compensating term originating from the expansion of 
^(n). 

4 Correction to P states 



For P states, Fg(nP,e) assumes the functional form 

9 V n 



F 6 (nP,e) = - (1--2 ) lne + iV(nP). (16) 
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We obtain the results, 



N(2P) = 0.003 300 635(1), N(SP) = 0.003 572 084(1) , 

AT(4P) = -0.000 394332(1), JV(5P) = -0.004 303 806(1) 

JV(6P) = -0.007496 998(1), JV(7P) = -0.010014614(1) 

jV(8P) = -0.011999 223(1). 



;i7) 



The result for iV(2P) was previously obtained in [9, Eq. (4.159)]. Specifically, the contribution 
due to this term entered into the result for F$h given in [6, Table I]. Here, we present the 
generalization of the result to higher principal quantum numbers. Observe that the correction 
changes its sign as n is increased. In analogy to the S states, we use a least-squares fit with a 
functional form 

JV(nP) « a P + bp/n + c P /n 2 , (18) 
which results in the fit-parameter values 

ap = -0.030 , bp = 0.170 , c P = -0.206 . (19) 

The excellent agreement between the fit and the numerical data is represented in Fig. 2. 

N(nP) 

0.008 
0.0 
-0.008 
-0.016 



2 3 4 5 6 7 8 



Figure 2: The analogue of Fig. 1 for P states. 



The spin-dependent high-energy part for P states is 

A^ H) (nP„ e) = ^ (Zaf ^ H 6 (nP jt e) . (20) 

Inspired by the effective treatment of radiative corrections based on form factors [10, Ch. 7], we 
take into account the contribution due to hard virtual photons by the replacement 

5V ^ 5V F 1 (-q 2 ) , (21) 

but this does not give the complete result. We also have to consider the perturbation of the 
spin-dependent interaction 

J_ (i7 . § E) - -L (i 7 • SE) F 2 (-q 2 ) , (22) 
2m 2m 
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where 5E is the electric field generated by the perturbating potential SV, i.e. 

ir(Za) 



5E 



-lq 



(23) 



in momentum space, and 7 is a three- vector whose elements are the spatial Dirac matrices [10]. 
Both of the above replacement prescriptions (21) and (22) are dictated by the modified Dirac 
Hamiltonian as given in [2, Eq. (3)]. 

We first consider the correction due to the charge form factor F\ as given in Eq. (21). Because 
we consider a priori only the one- loop effect, we may take F± in the one- loop approximation and 
expand only up to the order of q 2 . A suitably chosen perturbating potential may be used for 
an effective treatment of further loops, as discussed in Sec. 5 below. Formulas relevant to the 
electron charge form factor F\ can be found in Eqs. (5), (10), and (31) of [2]. We find 



a 



Za 
3^? 



In 



m 
27 



+ 



11 

24 



A(|Vnp(r)| : 



r=0 



(24) 



where ip n p( r ) is the nonrelativistic wave function of the nP state, and the matrix element reads 



A(|^nP(r)| 2 ) 



r=0 



2 (Zafm 5 

3 ~ 



7r n° 



1 



Therefore, the result for the scaled high-energy part due to the F\ form factor reads 
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■n. 



m 



2c 



11 



24 



(25) 



(26) 



Now we turn to the spin-dependent correction proportional to the magnetic form factor F2 
as given in Eq. (22). Again, we may employ the one-loop approximation and take F2 at zero 
momentum, where its well-known Schwinger value reads i*2(0) = a/(2n). We find 



A£jf' 2) (nPi, e) = F 2 (0) (nP,- | 7 ■ q\ nP-) , 

where the expectation value of the momentum space operator 7 • q reads 



7T (Za) 



(nPj I7 • q\ nPj) = i 



n 2 -l 



r=0 



irn 



(Za) m 5 jA/2 , 



(27) 



(28) 



i.e. the matrix element vanishes for P3/2 states [please observe the missing factor ir in the 
denominator of the right-hand side of Eq. (29) of [2]]. We thus obtain 



Hf\nV^e) 



1 



1 



n- 



(29) 



in analogy to the result obtained in Eq. (30) of Ref. [2] . 

Adding the low-energy part (16), and the two high-energy contributions (26) and (29) due to 
F\ and F2, respectively, we obtain 



AEsinPj) = AEl L) (nP,e) + £ AE^(nP j} e) = ^ (Za) 6 ^ F & (nP j 



(30) 



fc=l,2 



where the function Fg(nP) is independent of e and has the form 
F^nPj) = F e (nP,e)+ ^ H^^nP^e) 



fc=l,2 



(" , ((^)- ln(2) + i + i^' 2 ) + Ar(nP) - 



(31) 



The numerical values of N(nP) are given in Eq. (17). 
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(d) (c) (f) 



Figure 3: Some of two-loop diagrams involving the bound-state elec- 
tron self-energy which contribute to the logarithmic coefficient Bq\. All 
of these (and others) were taken into account in Ref. [3]. The double 
line in the diagrams denotes the bound electron (Dirac-Coulomb prop- 
agator). Figs, (a)-(c) contribute also to the triple logarithm of order 
or (Zaf ln 3 [(Za)~ 2 ] m, i.e. to the -E>63( n S)-coefficient. The double loga- 
rithm of order a 2 (Za) 6 ln 2 [(Za)~ 2 ] m originates from diagrams (a)-(e). 
General formulas valid for B^inS) and Bq3(hS) of arbitrary n are given 
in Eqs. (48), (49) and (51) of Ref. [3]. The single logarithm 5 61 (nS) 
involves many more contributions, for example those originating from 
the diagram (f). This Feynman diagram is considered in Eq. (40) of [3], 
together with contributions from other diagrams involving vacuum po- 
larization effects. Complete results for Bqi(uS) are given in Eq. (32). 



5 Conclusions 

We summarize the results of this Letter. 

• Dependence of the Bqi coefficient on the principal quantum number. Currently, 
one of the most important limiting factors to a further progress of the theory of the S 
state Lamb shift is the understanding of higher-order binding corrections to the two-loop 
bound-state self-energy [11]. The theory of S states is important for the deduction of the 
Rydberg constant, and it is also necessary to investigate higher excited states because the 
frequencies of more than one transition have to be theoretically known in order to infer 
the fundamental constants (see [11] and references therein). 

Our result for iV(lS) in Eq. (12) differs slightly from previously published results [see 
Eq. (21) of [3]]. General formulas, whose structure is valid for Bqi(uS) of arbitrary n have 
been given in Ref. [3, Eqs. (50) and (52)]; these are relevant to the sum of two-loop self- 
energy and vacuum-polarization effects. Some of the contributing Feynman diagrams are 
shown in Fig. 3. However, the quantity N(nS) which enters into the expression for Bqi(uS) 
has been known only for n = 1. Here, we take into account the slightly shifted valued of 
iV(lS) as well as the results presented in Eq. (12) for n > 1. We finally obtain the following 
results for Bei(nS) defined in Eq. (1), 

S 6 i(lS) = 50.344005(1), B 6 i(2S) = 42.447669(1) , 
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J3 6 i(3S) = 40.289637(1), B 61 (4S) = 39.294929(1) , 



B 6 i(5S) = 38.722048(1), B fl i(6S) = 38.348 805(1) , 

B 61 (7S) = 38.085860(1), B 6 l(8S) = 37.890 354(1) . (32) 

This completes the calculation of logarithmic two-loop corrections to S states of order 
a 2 {Zaf W[(Za)- 2 ]m (j = 1,2,3). 

• Clarification of a specific intermediate step in the calculation of the muonium 
hyperfine splitting. In lowest order, the hyperfine splitting Hamiltonian for an electron 
bound to an atomic nucleus with magnetic moment /x n is given by [10, p. 79]: 

Vhi s (r) = cr c • [V x ( Mn x V] 



2m 4vrr 

2 



<irm 



»n " {H n ■ V) V 



r 



-> _£_(cr e .„J *»)(,.), (33) 
3 m 

which is a delta-like potential as in (3). The angular averaging V$ Vj — > (1/3) 5ij V 2 [see 
the transition from the second to the third line of Eq. (33)] is valid when Vhf s (r) is evaluated 
on S states. Therefore, the quantity N(nS) appears naturally in the evaluation of radiative 
corrections to the hyperfine splitting of S states, for example in muonium [7, 12]. 

The results for N(1S) given in Eq. (12) differs slightly from the corresponding previously 
published value for the low-energy part of the muonium hyperfine splitting as given in [7]. 
The difference is 



2/31 vr 2 
3 V 36 ~ 12 



0.0257627... (34) 



While the new result reported here does not affect the final result for the hyperfine splitting, 
it explains the discrepancies between the intermediate results of [7] and those of Ref. [12]. 

Vacuum-polarization induced correction to the self-energy. Our results for N(nP) 
in Eq. (17) and the total results in Eqs. (30) and (31) for the energy shift due to the radia- 
tive correction to a delta-like potential can be used in order to evaluate the contribution 
to the energy of a P state due to the diagram in Fig. 3 (d). This is a combined "self- 
energy vacuum-polarization" correction for P states generated by a vacuum polarization 
correction to the Coulomb exchange in the bound electron propagator within the one-loop 
self-energy diagram. The lowest-order one-loop vacuum-polarization potential reads [10, 
p. 327] 

4 a K{Za)5^(r) 

VV.P. (r ) = - — 2 ' ( 35 ) 

15 7r m z 

which is proportional to 8V defined in Eq. (3). Our results in Eqs. (17), (30) and (31) 
contribute to the two-loop coefficients Bqi and Bqq for P states as defined in [2]. 

Nuclear finite-size correction to the self-energy. For an atomic nucleus whose root- 
mean-square charge radius is small compared to the Bohr radius, the effect of the nuclear 
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finite size can be incorporated as a form-factor correction to the Coulomb interaction in 
full analogy to Eqs. (21) and (22). In this case the appropriate form factor reads 



F(-q 2 ) = 1 - q 2 ^ (36) 

where (r 2 ) is the mean square radius of the charge distribution of the atomic nucleus. The 
Dirac delta-like finite-size "potential" therefore reads 



^ = >^*JZ4pa, (37) 



which is again proportional to 8V as defined in Eq. (3). Therefore, the results in this 
Letter can be used for an evaluation of the finite-size correction to the self-energy in 
atomic systems with a low nuclear charge number. Our result in Eq. (17) for the case 
n = 2 has been obtained previously in [9]. For the 2Pi/ 2 and 2P 3/ / 2 states, our calculations 
are also in agreement with Eqs. (16) and (17) of Ref. [13]. Our results in Eqs. (17) and 
(31) above represent generalizations of this previous work to higher excited states. 

Work on the nonlogarithmic term Bqq for S states is currently in progress and will be presented 
elsewhere, together with improved theoretical values for the Lamb shift of S states as derived 
from the results presented in this Letter. 
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